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We present a systematic assessment of the structural properties, the electronic density of states, 
the charge densities, and the phase stabilities of AgInSe 2 and AuInSe 2 using screened exchange 
hybrid density functional theory, and compare their properties to those of CuInSe 2 . For AgInSe 2 , 
hybrid density functional theory properly captures several experimentally measured properties, in¬ 
cluding the increase in the band gap and the change in the direction of the lattice distortion pa¬ 
rameter u in comparison to CuInSe 2 - While the electronic properties of AuInSe 2 have not yet been 
experimentally characterized, we predict it to be a small gap (« 0.15 eV) semiconductor. We also 
present the phase stability of AgInSe 2 and AuInSe 2 according to screened-exchange density func¬ 
tional theory, and compare the results to predictions from conventional density functional theory, 
results tabulated from several online materials data repositories, and experiment (when available). 
In comparison to conventional density functional theory, the hybrid functional predicts phase sta¬ 
bilities of AgInSe 2 in better agreement with experiment: discrepancies in the calculated formation 
enthalpies are reduced by approximately a factor of three, from « 0.20 eV/atom to ~ 0.07 eV/atom, 
similar to the improvement observed for CuInSe 2 . We further predict that AuInSe 2 is not a stable 
phase, and can only be present under non-equilibrium conditions. 


I. INTRODUCTION 

Chalcopyrite materials are important for a wide range 
of optoelectronic devices, for instance CuInSe 2 is well- 
known for its use as an absorber layer for applications in 
solar energy conversion. Thin-film photovoltaics based 
on the compound Cu(In,Ga)Se 2 (CIGS) have demon¬ 
strated the highest efficiencies to date among thin film 
polycrystalline materials 1 . In addition to CuInSe 2 , other 
chalcopyrite materials such as AgInSe 2 have also recently 
become of interest 2 . Silver, which has the same num¬ 
ber of valence electrons as copper, can be alloyed into 
CIGS in order to increase the band gap without causing 
structural disorder 3,4 . Although less studied to date, the 
chalcopyrite AuInSe 2 has been suggested as a candidate 
material to observe a three-dimensional topologically in¬ 
sulating state 5,6 . In addition, AuInSe 2 is reported to 
form during the vapor-liquid-sold growth of indium se- 
lenide nanowires from gold catalyst particles as well as 
possibly at the interface between gold/indium selenide 
metal-semiconductor junctions 7 9 . 

However, while the properties of the chalcopyrite 
CuInSe 2 have received a lot of attention in the 
literature 10-12 , related materials such as AgInSe 2 13,14 
and AuInSe 2 5 have received comparatively less. These 
and other emerging materials may also be of interest for 
a variety of related applications, but their properties are 
less well understood. It is unknown how effective alloy¬ 
ing AgInSe 2 into CuInSe 2 will be as a means to modify 
the band gap for photovoltaic applications and how the 
presence of Ag may affect transport or defect proper¬ 


ties. Also, it is of interest to know what the properties 
of AuInSe 2 would be, if it were possible to synthesize the 
compound in the laboratory. Since there are no reports 
of the targeted synthesis of this compound, guidance on 
how it may be possible to synthesize the material may 
be useful. Therefore, in this work we present a com¬ 
prehensive analysis of the electronic structure and phase 
stability of the selenide chalcopyrite materials AgInSe 2 
and AuInSe 2 using screened exchange density functional 
theory. 

It is known from the literature that obtaining an 
accurate first-principles description of the electronic 
structure of CuInSe 2 is challenging, due to the mixed 
ionic/covalent nature of the bonding and the need to 
properly describe the p-d hybridization in the solid 15,16 . 
For instance, within the local density approximation or 
the generalized gradient approximation in density func¬ 
tional theory (DFT), the electronic band gap of CuInSe 2 
(1.04 eV in experiment) closes. Hybrid density function¬ 
als have been shown to improve the description of the 
electronic structure and phase stability of CuInSe 2 10,11 , 
obtaining a better description of the electronic structure 
at a reasonable computational cost. 

Based on these considerations, we use hybrid screened- 
exchange DFT 17 19 and verify that it performs well on 
AgInSe 2 in addition to CuInSe 2 for electronic and ther¬ 
modynamic properties. The hybrid functional HSE06 
correctly captures the increase in the band gap of 
AgInSe 2 compared to CuInSe 2 , which is consistent with 
the predicted description of p-d hybridization and struc¬ 
tural parameters. We further demonstrate that the phase 
stability of AgInSe 2 , calculated with the hybrid func- 
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tional, compares well to experiment and offers similar 
improvements over conventional DFT based on PBE 20 as 
observed for CuInSe 2 - Having verified the performance 
of the hybrid for the known chalcopyrites, we use it to 
predict the electronic and thermodynamic properties of 
AuInSe 2 , which are as of yet unreported. We predict 
AuInSe 2 to be a small gap semiconductor (in contrast to 
the predictions from data mining approaches 21 ), which 
can only be present as a non-equilibrium phase. This 
suggests that synthesis of AuInSe 2 likely needs to pro¬ 
ceed through non-equilibrium or potentially high pres¬ 
sure routes. 

In Section II, we describe our methodology in detail. In 
Section III, we present the electronic structure of AgInSe 2 
and AuInSe 2 in comparison to CuInSe 2 . Section IV pro¬ 
vides a discussion of phase stabilities of the compounds. 


II. METHODOLOGY 


The DFT 22,23 calculations presented here were car¬ 
ried out within the Vienna A b Initio Software Package 
(VASP) 24,25 , using both the PBE 20 and the HSE06 19 ap¬ 
proximations to the exchange-correlation potential. We 
used projector augmented wave (PAW) pseudopoten¬ 
tials 26,27 to represent the core electrons; we chose Ar, Kr, 
and Xe cores respectively for Cu, Ag, and Au, which keep 
the outer d electrons in the valence. For each compound, 
the Kohn-Sham orbitals are expanded in a plane wave 
basis set with sufficient energy cutoff and k -point sam¬ 
pling of the Brillouin zone so that all computed param¬ 
eters are converged to the number of significant figures 
shown. Total energies are reported for each compound 
at the DFT-optimized lattice constants according to each 
functional; when applicable internal atomic coordinates 
are relaxed so that all forces on each ion are < 0.01 eV/A 
for PBE as well as for HSE06. 

Hybrid functionals are a class of approximations to the 
unknown exchange-correlation energy functional that in¬ 
corporate a portion a of exchange from Hartree-Fock the¬ 
ory. The HSE06 functional 19 used here uses a screened 
Coulomb potential to calculate the exchange. Although 
there have been attempts to match the band gap to the 
experimental value by adjusting the range separation pa¬ 
rameter uj and/or exchange mixing parameter <a, in this 
work we set the simulation parameters to the suggested 
values a = 25%, uo = 0.2 A -1 for benchmarking. In 
the following, in addition to AgInSe 2 and AuInSe 2 , for 
completeness and to facilitate comparison we also in¬ 
clude results for CuInSe 2 (which have been previously 
reported 10,11,28 ); our results for CuInSe 2 are in good 
agreement with others. 

The visualization program VESTA is used to an¬ 
alyze and visualize the computed electron density 
distribution 29 . 


III. STRUCTURAL PARAMETERS AND 
ELECTRONIC STRUCTURE 


A. Structural Parameters 

The chalcopyrite structure ABC 2 (space group I42d) is 
closely related to the zinc blende lattice typical of binary 
semiconductors. However, two cation species (A, B ) are 
present in an ordered arrangement on the usual cation 
sublattice. As a result, each anion C has two A and two 
B cations as nearest neighbors, which leads to a sym¬ 
metry breaking tetragonal distortion. The conventional 
unit cell of the chalcopyrites is thus double the conven¬ 
tional unit cell of the zinc blende lattice, with two of the 
latter stacked together and c % 2a. Additionally, the 
anion C typically assumes an equilibrium position closer 
to one pair of cations than to the other (bond lengths 
are unequal). Therefore, the XInSe 2 (X = Cu, Ag, Au) 
chalcopyrites can be characterized by three structural pa¬ 
rameters: the lattice constants a and c and the anion dis¬ 
placement u. The u parameter represents how the anion, 
Se, is displaced from its ideal tetrahedral site: 

_ 1 1 d? x _ Se - d 2 In _ Se 

n — “h 0 

4 a 2 

where dx-Se and di n -Se are the bond length of X-Se and 
In-Se, respectively. The value u = 0.25 for tetrahedral 
bonding, with u < 0.25 denoting displacement towards 
the X atom. Since u is related to the degree of orbital 
hybridization between different atoms, many calculated 
properties such as band gap are sensitive to it 12 ’ 15 46 > 

The structural parameters for XInSe 2 , (X = 
Cu,Ag,Au) calculated by DFT using both the PBE and 
HSE06 and the corresponding experimental values are 
summarized in Table I. For X = Cu and Ag, the lattice 
constants a and c estimated by PBE are larger than ex¬ 
perimental data, showing deviations of around 2% from 
the experimental value. HSE06 improves a and c, bring¬ 
ing the values closer to experiment (within 1%). For 
CuInSe 2 , our calculated PBE and HSE06 parameters are 
in agreement with previously reported values 10,11 . For 
AgInSe 2 , our PBE results are consistent with previously 
reported values 13 . 

For AgInSe 2 the description of the distortion u is also 
improved in HSE06, similar to the improvement observed 
for CuInSe 2 - According to Table I, the deviation of u 
from experiment is large in PBE (^ 3%) but HSE06 im¬ 
proves it to within 0.3%. Notably HSE06 properly cap¬ 
tures the trend u < 0.25 (Se atom shifts towards Cu) for 
CuInSe 2 30 but u > 0.25 (Se atom shifts away from Ag) 
in AgInSe 2 31 . The difference in the direction of the dis¬ 
tortion may be due to the smaller size of Cu atoms com¬ 
pared to Ag atoms and/or a stronger interaction matrix 
element betweeen the Cu-s/Se-p core states in compari¬ 
son to Ag-s/Se-p. Since X-Se spacing is related closely to 
the degree of hybridization between the transition metal 
d orbitals and the 4p orbitals of the chalcogen, it is en¬ 
couraging that the match to experiment improves. 
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Table I. Calculated structural parameters a, c/a and u , and band gaps of VInSe 2 (X — Cu, Ag, Au). Percent differences 
compared to experimental data (when available) are shown in parentheses. 




a (A) 



c/a (A) 



u 



Band gap (eV) 


PBE 

HSE06 

Exp. 

PBE 

HSE06 

Exp. 

PBE 

HSE06 

Exp. 

PBE 

HSE06 

Exp. 

CuInSe 2 

5.88 

(1.8) 

5.83 

(0.8) 

5.78 a 

2.01 

(-0.1) 

2.01 

(-0.2) 

2.01 a 

0.22 

(-3.7) 

0.23 

(0.3) 

0.23 a 

0 

0.83 

(-19.8) 

1.04 c 

AgInSe 2 

6.20 

(1.6) 

6.16 

(1.0) 

6.10 b 

1.95 

(1.6) 

1.93 

(0.8) 

1.92 b 

0.25 

(-2.4) 

0.26 

(0.3) 

0.26 b 

0 

0.97 

(-22.1) 

1.24 c 

AuInSe 2 

6.14 

6.07 

- 

1.99 

1.98 

- 

0.24 

0.25 

- 

0 

0.16 

- 


a Reference 30 
b Reference 31 
c Reference 32 
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Figure 1. Electronic band structure of CuInSe 2 according to 
DFT-PBE (left) and DFT-HSE06 (right), in comparison to 
recent ARPES measurements (middle). Adapted with per¬ 
mission from Ref. [33]. Copyrighted by the American Physi¬ 
cal Society. 



Table I also provides the calculated structural param¬ 
eters for AuInSe 2 according to PBE and HSE06. While 
we cannot compare it to experimental data, it would be 
of interest to know whether as in the cases of V=Cu, Ag 
the use of the hybrid improves the description, or whether 
for X=Au the system is sufficiently delocalized that the 
PBE description is more appropriate. Based on the cal¬ 
culated parameters, we predict that the lattice constants 
of AuInSe 2 lie in between those of CuInSe 2 and AgInSe 2 , 
and that u ~ 0.25, suggesting less symmetry breaking 
and a more tetrahedral structure. 


B. Electronic Structure and Density of States 


In comparison to the corresponding binary II-VI semi¬ 
conductors, the chalcopyrites exhibit a suppressed band 
gap: for instance in CuInSe 2 the gap is E g = 1.04 eV 32 , 
much lower than its II-VI counterpart ZnSe for which 
Eg = 2.68 eV 32 . The reason for the suppressed gap can 
be understood from Fig. 1, which shows the structure of 
CuInSe 2 valence bands according to PBE and HSE06, in 
comparison to recent angle-resolved photoemission spec¬ 
troscopy (ARPES) measurements 33 . The differences be¬ 
tween the PBE and HSE06 band structure of CuInSe 2 
have been noted previously 28 ; we summarize the key con¬ 
cepts here to observe similarities and differences with 
AgInSe 2 and AuInSe 2 - The projected density of states 
(PDOS) of the CuInSe 2 valence bands onto the atomic 
states are also shown in Fig. 2 (leftmost column), which 
show the nature of the electronic structure as (1) upper¬ 
most VBs: antibonding Cu 3d/Se 4 p states, (2) 2 eV 
below VBM: a set of non-bonding Cu 3d states, and (3) 
4 eV below VBM: the bonding Cu 3d/Se 4p levels, and 
(4) 6 eV below VBM: bands derived from interactions be¬ 
tween In and Se atoms. According to Fig. 1, for CuInSe 2 , 
HSE06 improves the orbital spacing and bandwidths of 
the valence bands in comparison to PBE (the relative 
positions are better matched to experiment). The cal¬ 
culated band structures for AgInSe 2 and AuInSe 2 (not 
shown) are qualitatively similar to that of CuInSe 2 in 
Fig. 1 in terms of the orbital ordering, but the orbital 
spacing and bandwidths vary. 

For CuInSe 2 , interactions between the non-bonding Cu 
3d states and the antibonding Cu 3d/Se 4 p levels towards 
the top of the valence bands in Fig. 1 leads to a repul¬ 
sive p-d interaction that pushes the antibonding states 
upwards in energy, which is largely responsible for the re¬ 
duced band gap 15,16 . Capturing this interaction properly 
is necessary to accurately position the bands and obtain a 
good estimate of the band gap. For instance, for CuInSe 2 
(Fig. 1) according to PBE the 3d non-bonding levels are 
too high in energy, resulting in an “overhybridization” 
with the upper valence bands. The overhybridization 
pushes the upper VBs further upwards and is largely re¬ 
sponsible for the closure (in fact, the inversion 34 ) of the 
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Figure 2. Density of states (DOS) for CuInSe 2 (left), AgInSe 2 (middle), AuInSe 2 (right) calculated by DFT-PBE (upper) and 
DFT-HSE06 (lower). The black lines are total DOS, and red, blue and green lines are projected DOS of Cu/Ag/Au, In, and 
Se respectively. The Fermi level is shown by the dashed grey line. 


gap in PBE, see Fig. 2. The use of HSE06 drops the 
position of Cu 3d non-bonding orbitals and opens the 
gap to E^ se06 = 0.83 eV, in better agreement with the 
experimental value of 1.04 eV 32 . 

To see what changes in the electronic structure of 
XInSe 2 for X=Ag, Au in comparison to X=Cu, the to¬ 
tal and projected density of states (DOS) of the XInSe 2 
compounds are shown in Fig. 2. The dotted line indicates 
the top of the valence band in all cases. The black lines 
are the total DOS, and red, blue, and green indicate the 
transition metal X, In, and Se states respectively. For 
AgInSe 2 the first observation is that the HSE06 band gap 
is larger than the PBE gap (as expected). In PBE the 
gap is closed (in fact, inverted), while in HSE06 it opens 
to 0.97 eV. The opening of the gap in AgInSe 2 in HSE06 
occurs for similar reasons as discussed above for CuInSe 2 : 
the DOS in Fig. 2 (middle column) shows that the non¬ 
bonding 4 d levels are shifted downwards in HSE06, re¬ 
ducing the p - d repulsion, and opening the gap. In 
addition, according to HSE06 the gaps are E g = 0.83 
eV for CuInSe 2 and E g = 0.97 eV for AgInSe 2 , in rea¬ 
sonable agreement with the experimental values (1.04 eV 
and 1.24 eV respectively). Therefore, HSE06 is also able 
to capture the increase in the gap when Cu is replaced 
with Ag (discussed further later). 

The encouraging performance of the hybrid functional 
for AgInSe 2 as well as CuInSe 2 suggests that it may 
be well-suited for AuInSe 2 as well. Interestingly we 
find that the trend described above persists somewhat 
for AuInSe 2 : the gap is closed (inverted) in PBE, but 
E g ~ 0.16 eV in HSE06. Thus, although the gap opens 
in HSE06, it does so to a lesser extent than the other 
chalcopyrites. To understand why, from Fig. 2 as X 
changes from Cu -T Ag Au, the relevant d orbitals 
become more delocalized in nature and exhibit a larger 


bandwidth in the DOS. Accordingly the degree of mixing 
between the non-bonding X-d and antibonding X-dj Se 
4 p states within the upper valence bands increases. The 
more delocalized nature of the 5 d bands is already bet¬ 
ter described in PBE, compared to the more localized 
3d and 4 d bands. Correspondingly, when HSE06 is used 
in AuInSe 2 , there is no need for as large a “correction” 
to the position of the nonbonding Au 5 d states and the 
HSE06 and PBE DOS (and thus, the gap) for AuInSe 2 
are more similar. 

Overall the HSE06 gap of AgInSe 2 (0.97 eV) is larger 
than that of CuInSe 2 (0.83 eV), and then quickly drops 
to 0.16 eV for AuInSe 2 - The increase in the gap as Cu -T 
Ag is consistent with experiment, but the non-monotonic 
trend as Cu —> Ag —>> Au is different from typical expec¬ 
tations of a reduced gap as isovalent elements in a com¬ 
pound are successively exchanged by descending along a 
column of the periodic table 35 . These trends can be bet¬ 
ter understood from the calculated distortion parameter 
u in Table I and the HSE06 charge densities shown in 
Fig. 3, plotted on a plane that contains X, In, and Se 
atoms. The scale bar and the color map isosurface levels 
are the same for all three plots, for direct comparison of 
the charge distribution in the three materials. 

According to Fig. 3, both ionic and covalent bonds are 
present along the In-Se and X-Se bonds, respectively, in 
the XInSe 2 system. Electrons are drawn towards Se and 
away from In in the In-Se bond, reflecting a more ionic 
nature, while the more evenly distributed charge in the 
X-Se bond indicates a more covalent nature. For the X- 
Se bond, arrows indicating the width corresponding to 
a particular isosurface value for all three materials are 
marked in Fig. 3. The arrow is widest for AuInSe 2 , nar¬ 
rowest for AgInSe 2 , and intermediate for CuInSe 2 - This 
ordering is consistent with the relative band gaps in the 
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Figure 3. Electronic charge density distribution of CuInSe 2 (left), AgInSe 2 (middle), AuInSe 2 (right), across a plane that 
contains Cu/Ag/Au, In and Se atoms according to DFT-HSE06. The contour lines of isosurfaces are spaced logarithmically 
and the arrows are marked corresponding to the same isosurface value for all three materials. 


compounds, as well as the differences in the u parameter 
and varying degrees of p - d hybridization for all mate¬ 
rials in Fig. 2. As u < 0.25 in CuInSe 2 , the Cu and Se 
atoms are closer together, and orbital overlap and p-d re¬ 
pulsion is greater in this material, resulting in a smaller 
band gap. By contrast, for AgInSe 2 where u > 0.25, the 
Se atom is further away from Ag and closer to In. Less 
direct orbital overlap and p-d repulsion within the upper 
valence bands results in a larger gap for this material. 
For AuInSe 2 , u = 0.25 suggests a more typical tetrahe¬ 
dral bonding. The high degree of charge sharing in the 
AuInSe 2 results from the larger atomic size of Au and 
the 6d orbitals, and the more delocalized nature of these 
orbitals results in smaller p-d repulsion and generally 
more delocalized states and metallic like character in the 
solid. 


IV. PHASE STABILITIES AND ACCURACY 
ASSESSMENT 


Having characterized the electronic structure of the 
materials, we now turn to the formation enthalpies and 
phase stabilities. Not only is this analysis directly of rel¬ 
evance for the chalcopyrites themselves, but it also sug¬ 
gests the degree of accuracy attainable in the calcula¬ 
tion of phase stabilities in general and thus has implica¬ 
tions for the long-standing goal of computational materi¬ 
als prediction and discovery 36 . Therefore, in the follow¬ 
ing the phase stability according to PBE and HSE06 are 
compared to each other and experiment (when available). 
Results for CuInSe 2 phase stability have been previously 
reported 10,11 , we show them here again for comparison 
purposes. For completeness, we also compare our com¬ 
puted phase stabilities to those that are available in sev¬ 
eral online materials repositories 37-39 to assess patterns 
and trends in the results. 


A. Formation Enthalpies 

In a thermodynamic framework, phase stability dia¬ 
grams determine the regions of chemical potential phase 
space for which a particular compound is stable with re¬ 
spect to the formation of all possible elemental solids and 
other binary or ternary compounds. Different growth or 
environmental conditions are accommodated by the set 
of chemical potentials pi for each atomic species i that 
is present in the compound, by assuming that each is in 
equilibrium with a physical reservoir such as a gas or a 
bulk phase. At thermodynamic equilibrium the chemical 
potentials satisfy 



where the summation is carried out over all atomic 
species present in the compound, e.g. PcuinSe 2 = 
PCu + Hin + 2 pse- By referencing the chemical potential 
of each atomic species in the compound to that of the 
bulk elemental solid or gas phase so that pi = p° + Ap i: 
then 

l^comp = ^ ^ l^i T ^ ^ • (3) 

i i 

The formation enthalpy of the compound, 

Aff = ^A Mi , (4) 

i 

is the change of enthalpy when the compound is formed 
from its constitutent elemental phases. Negative AH 
denotes stable compounds. Note that we neglect TS and 
PV contributions to the enthalpy here. 

The calculated formation enthalpies of the chalcopy- 
rite materials according to PBE and HSE06 are shown 
in Table II, in comparison to the experimental values. For 
AgInSe 2 , the enthalpy of formation from PBE is under¬ 
estimated by more than 30% of the experimental values 
(similar to the underestimate observed for CuInSe 2 ). The 
underestimate occurs because the elemental compounds 
Cu or Ag, In, and Se are better described in PBE than 
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the chalcopyrites themselves. Since the chalcopyrites and 
the elemental solids exhibit different degrees of metallic, 
covalent, or ionic character, PBE suffers from an incom¬ 
plete cancellation of errors. The incomplete cancellation 
of errors arises because PBE is not able to describe ma¬ 
terials with different degrees of electron localization on 
equal footing 38,41,42 . 

According to Table II, HSE06 also underestimates the 
formation enthalpy of AgInSe 2 , but now the error is re¬ 
duced to « 11%. Therefore, it appears that HSE06 im¬ 
proves the description overall (and the level of improve¬ 
ment is quite similar to that observed for CuInSe 2 ). For 
comparison, we also show the formation energies accord¬ 
ing to NREL Materials Database 38 , the Open Quantum 
Materials Database (OQMD) 39 , and Materials Project 37 . 
The results from these online repositories mostly match 
well our PBE calculations, as expected since similar cal¬ 
culation protocols are used. For example, the Materials 
Project uses DFT to evaluate the total energy of com¬ 
pounds and mixes GGA and GGA+U for the exchange- 
correlation functional 41 for accurate prediction of forma¬ 
tion enthalpies. The NREL Materials Database forma¬ 
tion enthalpy is different from the other two, as a result 
of the Fitted Elemental Reference Energies (FERE 38 ) ap¬ 
proach. Interestingly in this case this brings the forma¬ 
tion enthalpies more into accord with our hybrid DFT 
results and closer to the experimental values. 

For AuInSe 2 , we predict formation enthalpies of -1.23 
eV/fu in PBE and -1.59 eV/fu in HSE06; results for this 
material are not available in the repositories or experi¬ 
ment. This suggests that the compound is favorable to 
form from isolated elemental solids Au, In, and Se. How¬ 
ever, the formation enthalpy is lower than the predicted 
value for the other chalcopyrite materials. For this com¬ 
pound, we note that the PBE and HSE06 description 
are more similar to each other on an absolute scale, in 
comparison to the other two chalcopyrites. 


B. Formation Enthalpies of Competing Binary 
Phases 

A detailed comparison of formation enthalpies of the 
chalcopyrites and several competing candidate binary 
phases is shown in in Fig. 4. Fig. 4a shows absolute 


Table II. Calculated enthalpy of formation of XInSe 2 (X = 
Cu, Ag, Au) per formula unit. Percent differences compared 
to experimental data are shown in parentheses. 



PBE 

HSE06 Mat Proj 37 

NREL 38 

OQMD 39 Exp 

liTn Qn 

-1.78 

-2.40 

-1.77 

-2.32 

-1.80 -2.77 a 

Uumbe2 

(-35.6) 

(-13.4) 

(-36.1) 

(-16.3) 

(-35.0) 

AgInSe 2 

-1.71 

-2.22 

-1.72 

-2.08 

-1.73 -2.51 a 

(-31.8) 

(-11.4) 

(-31.5) 

(-17.1) 

(-31.1) 

AuInSe2 

-1.23 

-1.59 

- 

- 

- 


a Reference 40 


values of formation enthalpies, and Fig. 4b shows the ab¬ 
solute deviations from the experimental values, both on 
a per atom basis. The compounds shown in Fig. 4 were 
selected by using the known experimental formation en¬ 
thalpies as a guide to identify the most competitive alter¬ 
nate phases. For CuInSe 2 , Cu-Se binaries such as Cu 2 Se, 
CusSe 2 , CuSe, CuSe 2 were selected, and for AgInSe 2 , 
Ag 2 Se was included. The In-Se binaries, In 2 Se 3 , InSe, 
In 2 Se were incorporated in the phase space for both cases. 
Some of these phases are metallic (I^Se and all of the 
copper selenides except for Cu 2 Se) and others semicon¬ 
ducting (the chalcopyrites, Cu 2 Se, Ag 2 Se, AuSe, In 2 Se 3 , 
and InSe), see the DOS according to HSE06 presented in 
the Appendix. In some cases in Fig. 4 (AuInSe 2 , AuSe, 
and In 2 Se) the experimental enthalpies of formation are 
not included because we are unable to find reported val¬ 
ues. 

Several trends are evident from Fig. 4. First, as ex¬ 
pected, it is clear that PBE underbinds the compounds 
with respect to the solid elemental phases. On a per 
atom basis, the worst cases are the chalcopyrites CuInSe 2 
and AgInSe 2 themselves and the semiconductor Cu 2 Se 
for which the errors are all larger than 0.20 eV/atom. In 
fact, Cu 2 Se is predicted to be slightly unstable in PBE. 
Based on their density of states (see Appendix), these 
worst cases correspond to semiconductors for which the 
valence bands have a more substantial d orbital char¬ 
acter. For these materials, HSE06 brings the formation 
enthalpy closer to the experimental value although it still 
underbinds (after all, ultimately, HSE06 is still not self¬ 
interaction free). For the chalcopyrites, the degree of 
underbinding is reduced to 0.07-0.09 eV/atom while for 
Cu 2 Se the error is reduced to within 0.01 eV/atom in 
comparison to experiment. 

For the remaining semiconductors Ag 2 Se, In 2 Se 3 , and 
InSe where the valence bands have less d orbital char¬ 
acter, both PBE and HSE06 both do better. In these 
cases, the PBE results are closer to experiment to be¬ 
gin with, exhibiting deviations from experiment of 0.07 
eV/atom, 0.14 eV/atom, and 0.12 eV/atom respectively. 
Here again, HSE06 results improve the description fur¬ 
ther, bringing the deviation to within 0.03 eV/atom, 0.02 
eV/atom, and 0.01 eV/atom respectively. 

Considering the metallic compounds Cu 3 Se 2 , CuSe, 
and CuSe 2 , the PBE description is better as expected: 
the formation enthalpy is always within ~ 0.07 to 0.08 
eV/atom in all cases. In the case of Cu 3 Se 2 and CuSe, 
HSE06 improves the description to within 0.01 eV/atom, 
but it does not change the description of CuSe 2 - 

For comparison purposes, in Fig. 4 we also show the 
formation enthalpies according to Materials Project and 
the NREL database. As expected the Materials Project 
results agree well with our PBE results. It is interesting 
to observe that, with the exception of AuSe for which 
there is a large deviation, for most of these compounds 
the NREL Materials Database results, obtained using 
FERE 38 , match reasonably well the results of the hybrid 
calculations presented here. 
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Figure 4. Comparisons of formation enthalpies of the chalcopyrites and the binaries calculated by DFT-PBE and DFT- 
HSE06 with experimental data for ATnSe 2 40 , Cu-Se and In-Se binaries 43 , and Ag-Se binary 44 , (a) Absolute values of formation 
enthalpies in eV per atom, (b) Absolute deviations from the experimental values in eV/atom. 


C. Phase Stability Diagrams 

Finally, phase stability diagrams for CuInSe 2 , 
AgInSe 2 , and AuInSe 2 obtained from the chalcopyrite 
and competing phase formation enthalpies are shown in 
Figs. 5, 6, and Fig. 7. For CuInSe 2 and AgInSe 2 we com¬ 
pare PBE and HSE06 calculation results to experimental 
values directly; for AuInSe 2 the experimental values are 
not known so only calculation results are presented. Each 
point inside the triangles corresponds to a set of chemi¬ 
cal potentials satisfying thermodynamic equilibrium, e.g. 
A/icu + A/ij n + 2AyUse = A H Cu inSe 2 , and the constraint 
that Afa < 0 so that the compound is not unstable with 
respect to decomposition to elemental solids. This gives 
rise to the outer triangles in Figs. 5, 6, and 7, the size 
of which is governed by the value of AH for the chal¬ 
copyrites. The shaded region corresponds to the part 
of phase space for which the chalcopyrite is also sta¬ 
ble with respect to decomposition to competing binaries 
TA A Hi ^ A H})i nar y . 

The trends observed for AgInSe 2 in Fig. 6 are quite 
similar to those of CuInSe 2 in Fig. 5. Since from Ta¬ 
ble II, the formation enthalpy in PBE is small compared 
to experiment, the size of the PBE triangle is small but 
the HSE06 triangle is larger and closer to experiment. 
Boundaries marking where the compound is unstable 
with respect to a particular binary are represented as 
lines in Figs. 5 and 6, with blue lines representing X-Se 
binaries and red lines representing In-Se binaries. It is 
interesting that even in PBE, despite the smaller size 
of the triangles and the large degree of underbinding 


of AgInSe 2 , the compound is still predicted to be sta¬ 
ble in some portion of phase space (the same holds for 
CuInSe 2 ). We note that this may not always be the case: 
it is conceivable that if the degree of underbinding were 
slightly larger, the compound may be incorrectly pre¬ 
dicted to be unstable (a “false negative”). This possibil¬ 
ity highlights some of the challenges in stability predic¬ 
tion, since accurate phase stability diagrams rely on ac¬ 
curate calculation of formation enthalpies across a wide 
range of competing phases from metals to semiconduc¬ 
tors to insulators. The challenge is complicated by the 
fact that there are many materials that are almost, but 
not, thermodynamically stable. That, combined with the 
uncertainty in DFT makes the materials search problem 
very challenging. 

Finally, Fig. 7 shows the phase stability diagram for 
AuInSe 2 according to PBE and HSE06. Unlike the other 
two compounds, the line for In 2 Se 3 line is not shown 
because it is located below the triangle. Thus, accord¬ 
ing to both PBE and HSE06 AuInSe 2 is always unstable 
with respect to the decomposition into In 2 Se 3 and Au 2 Se, 
and we predict that under all circumstances it can only 
be present as a non-equilibrium phase. Given that both 
functionals predict this and that the In 2 Se 3 line lies so 
far below the triangle, we do not believe this is a case of a 
“false negative”, which would require that the error in the 
calculated formation enthalpies be much larger than ob¬ 
served for the other chalcopyrites. Therefore, synthesis of 
AuInSe 2 likely needs to proceed through non-equilibrium 
or potentially high pressure routes. 
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Figure 5. Phase stability diagrams of CuInSe 2 according to DFT-PBE (left), DFT-HSE06 (middle), and experimental data 
(right). The area below the blue lines or above the red lines indicate where CuInSe 2 is unstable with respect to formation of 
Cu-Se or In-Se binaries, respectively. The shaded area shows the chemical potential range for which CuInSe 2 is stable. 



Figure 6. Phase stability diagrams of AgInSe 2 according to DFT-PBE (left), DFT-HSE06 (middle), and experimental data 
(right). The area below the blue line or above the red lines indicate where AgInSe 2 is unstable with respect to formation of 
Ag-Se or In-Se binaries. The shaded area shows the chemical potential range for which AgInSe 2 is stable. 


V. CONCLUSIONS 

In conclusion, we have carried out (i) an in-depth, 
systematic comparison of the structural parameters 
and electronic structure of the chalcopyrites CuInSe 2 , 
AgInSe 2 , and AuInSe 2 using screened exchange hybrid 
density functional theory, and (ii) a comparison of the 
phase stability predictions in this material set accord¬ 
ing to PBE and HSE06 in order to assess the degree 
of accuracy attainable. Consistent with experiment, we 
find that the band gap of AgInSe 2 is larger than that of 
CuInSe 2 due to a reduced p-d repulsion in the valence 


bands. According to our HSE06 calculations, AuInSe 2 is 
predicted to have a small band gap ~ 0.16 eV. Regarding 
phase stability, for the known compounds CuInSe 2 and 
AgInSe 2 the calculated stabilities indicate that the com¬ 
pounds are stable in a range of chemical potential phase 
space in good agreement with experimental results. We 
further predict that AuInSe 2 is not stable, and can only 
be present as a non-equilibrium phase. A comparison of 
the phase stabilities suggests that for these chalcopyrites 
typical PBE errors in the calculation of formation en¬ 
thalpies can be as large as 0.20 eV/atom. However the 
use of HSE06 is observed to reduce this by approximately 
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Figure 7. Phase stability diagrams of AuInSe 2 from DFT- 
PBE (left), DFT-HSE06 (right). The area below the blue 
lines or above the red lines indicate where Au-Se or In-Se 
binaries are stable than AuInSe 2 . AuInSe 2 is not stable in 
either DFT functional for any value of chemical potentials. 


a factor of three. 
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Figure 8. Density of states (DOS) for binary compounds calculated by DFT-HSE06. The black lines are total DOS, and red, 
blue and green lines are projected DOS of Cu/Ag/Au, In, and Se respectively. The Fermi level is shown by the dashed grey 
line. 














































